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Abstract 

One-dimensional and quasi-one-dimensional strongly-coupled dusty plasma rings have been cre- 
ated experimentally. Longitudinal (acoustic) and transverse (optical) dispersion relations for the 
1-ring were measured and found to be in very good agreement with the theory for an unbounded 
straight chain of particles interacting through a Yukawa (i.e., screened Coulomb or Debye-Hiickel) 
potential. These rings provide a new system in which to study one-dimensional and quasi-one- 
dimensional physics. 
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One-dimensional (ID) chains of particles are an important model system for studying 
linear and nonlinear phenomena such as waves, instabilities, and diffusion [TH3] • A standard 
model is that of an unbounded straight chain of particles interacting through a given force. 
This system may represent a lumped model of a continuous system, such as a string, or may 
model the atomic lattice. Although infinite systems cannot be realized in practice, they can 
be approximated using large systems with periodic boundary conditions so that end effects 
are eliminated. One experimental system where it is possible to study strongly-coupled ID 
chains is dusty (complex) plasma. 

Dusty plasma is a system of microscopic material particles immersed in an electron-ion 
plasma. For typical laboratory conditions, the dust particles acquire a net charge q < 
and can be confined in a single-layer suspension above a horizontal electrode. Vertical 
confinement is due to the balance between the downward gravitational force and the average 
upward electrostatic force inside the plasma sheath. Horizontal confinement is often due to 
depressions in the plasma potential created by secondary electrodes. The interaction between 
dust particles is screened by the response of the electrons and ions so that the potential a 
distance r from a dust particle is given by a Yukawa (screened Coulomb or Debye-Hiickel) 
potential [U E] 



where Ad is the Debye screening length. In principle, the interaction length can be varied 
from long to short range allowing study of both "plasma" and "condensed matter" regimes 
[6]. Since the dominant friction force on dust particles is neutral gas (Epstein) drag, wave 
modes may be underdamped at low gas pressures. 

Previous experimental work on ID dusty plasmas has focused on straight chains confined 
in highly anisotropic two-dimensional (2D) biharmonic wells [THT2]. Such systems have 
several drawbacks. First, the lattice constant a is not a constant and increases at the ends 
of the chain. Second, the boundary conditions are not periodic, so that comparisons to 
unbounded theories are problematic. Third, it is difficult to make very long straight chains 
since ID chains are unstable against the zigzag when the lattice constant is below a critical 
value [12j . To overcome these drawbacks it may be possible to create a one-dimensional 
Coulomb ring [13] or dusty plasma ring [H] using a two-dimensional (2D) annular potential 
well. A dusty plasma ring will provide the most direct test of the theory for lattice waves in 
unbounded Yukawa chains [15J, as well as a system with which to explore phenomena such 
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as single-file diffusion [16] and the properties of quasi-one-dimensional systems |17j . 

In the present work we demonstrate that it is possible to create one-dimensional and 
quasi-one-dimensional dusty plasma rings in experiment, thereby providing a real system 
that models a chain of interacting particles with periodic boundary conditions. We measure 
the longitudinal and transverse dispersion relations for a 1-ring using thermal fluctuations, 
and confirm that the waves are well described by the theory of an unbounded ID Yukawa 
chain. 

We briefly review dusty plasma ring theory [Hj. The dusty plasma consists of n identical 
particles with mass m and charge q at positions {xi,yi} confined in a 2D annular potential 
well [13]. The potential energy of an isolated particle in this well is 

U we \\ (r) = ^muo (r - s) 2 , (2) 

where Uq is the single-particle radial oscillation frequency, r is a radial coordinate and s is 
the radius at which the well potential energy is a minimum. The total potential energy of 
a particle configuration is 
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where the second sum is the interaction potential energy for all unique pairs of particles and 
is the distance between particles % and j. The potential energy U is nondimensionalized 
by defining variables 

« = T— , o- = — , Q = — , Vi = — , Pij = ( 4 ) 
Ad r r r r 

where the characteristic length and energy are 
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to give 



In dimensionless form, time-independent configurations {£,i,r]i} are determined by three 
parameters: n, k and a. Here k is the Debye shielding parameter and a is the radius 
at which the confining well has its minimum. This model was shown [H] to have one- 
dimensional ring solutions (i.e., 1-rings). A 1-ring is unstable against the zigzag instability 



when the particle density exceeds a critical value, which may lead to the creation of quasi-lD 
systems consisting of concentric rings, e.g., 2-rings. 

For a ID chain with lattice parameter a the dimensionless lattice parameter is a = ajr^. 
The curvature of the ring is negligible when s>a and s ^> Xd- In this case, an appropriate 
shielding parameter is 

R = - — = «k. (7) 
Ad 

The dispersion relation for an undamped longitudinal (acoustic) wave in an unbounded 
straight ID Yukawa chain is [TS] 
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where is the longitudinal wave frequency and k is the real wavenumber. The dispersion 
for the transverse (optical) mode is [10] 

-± = 1 - _ £ _ (1 + [1 - cos (jfca)] , (9) 

where cut is the transverse wave frequency. The optical mode is a backward wave. These 
dispersion relations were previously found to be in good agreement with theoretical normal 
mode spectra for a dusty plasma 1-ring |14j . 

Our experimental setup is shown in Fig. [TJ Monodisperse microscopic (8.94 ± 0.09 fim 
diameter) melamine dust particles are levitated in an rf plasma discharge above the 
capacitively-coupled powered electrode in the DONUT (Dusty ONU experimenT) appa- 
ratus [7J H5H2T] . Dust confinement is due to a 2D annular potential well created by cutting 
an annular rectangular groove in the surface of the 89-mm powered electrode. The groove 
has an inner diameter of 15.9 mm, an outer diameter of 40.0 mm and a depth of 3.3 mm. 
The plasma sheath edge partially conforms to the groove, creating an annular depression 
in the equipotential surfaces above the groove. The dust particles are illuminated using a 
steady sheet of 635-nm laser light. Particle positions are recorded using a video camera. For 
these experiments 1025 frames of 1024 x 1024 pixels were recorded at 7.51 frames/s. 

Making a uniform dusty plasma ring [Fig. [2] requires the bottom of the potential well to 
be flat. This condition is more stringent than confinement in a parabolic well since a tilted 
parabolic well is still parabolic. To create a dust ring we adjusted the tilt of the electrode 
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Figure 1: Schematic of the experimental setup. An annular groove with rectangular cross section 
cut into the rf electrode produces an annular confining well for dust particles. When the particle 
spacing is above a critical value the particle configuration is one-dimensional with a ring topology. 

(actually the entire vacuum vessel) until the particles were distributed approximately uni- 
formly in a ring. However, at low pressures (< 20 mtorr) we observed an m = 2 azimuthal 
perturbation to the potential well. The exact characteristics of this perturbation depend on 
the rf power, suggesting it may be due to an rf standing- wave component. Such a pertur- 
bation cannot be eliminated by tilting the electrode. We used two grounded L-shaped rods 
inserted horizontally above the electrode to cancel this perturbation by varying the distance 
and orientation of the rods. 

An equilibrium configuration of a dusty plasma 1-ring is shown in Fig. [2]^a). The ring 
consists of n = 68 dust particles floating 9.2 mm above the top of the electrode. Discharge 
parameters were a neutral pressure of 18.9 mtorr argon, ~ 3 W forward rf power at 13.56 
MHz, and the dc self-bias on the electrode was -52 V. Wave modes are underdamped at 
this pressure. The camera resolution was 0.0317 mm/pixel . To determine the geometric 
center of the ring we fitted the particle coordinates to a circle j22] and then averaged over 
all frames. (Because of variable particle spacing the geometrical center and the center of 
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Figure 2: (a) Measured equilibrium positions of n = 68 particles confined in a two-dimensional 
annular well. The particles form a one-dimensional ring, (b) Distance a between nearest neighbors 
vs azimuthal angle. 

mass are not identical.) We found an average radius R = 11.73 mm, so that the average 
particle separation computed as the circumference over n is a = 1.084 mm. The dust ring 
lies closer to the inner edge of the groove. The Euclidean distance Tij between adjacent pairs 
of particles is shown in Fig. |2]^b). The average separation (ry) = 1.084 ± 0.008 mm agrees 
with the previous estimate of a. 

In Fig. [3]^a) we demonstrate that it is also possible using this confinement geometry to 
create a quasi-lD configuration of two rings (i.e., a 2-ring) with different radii, as predicted 
by theory [131 E]- m this case, we have a total of n = 87 particles, with 43 particles in the 
inner ring and 44 particles in the outer ring. The radius of the inner ring is 10.12 ± 0.01 
mm, and that of the outer ring is 11.24 ± 0.01 mm. For most of this two-ring the particles 
are in a zigzag configuration. However, near 0° the two rings "slip" past each other. The 
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Figure 3: (a) Quasi-one-dimensional configuration (a 2-ring) consisting of n = 87 particles arranged 
in two concentric rings with 43 particles in the inner ring and 44 in the outer ring, (b) Distance 
between particles in the inner ring (open circles) and in the outer ring (open diamonds) vs angle. 



distance between pairs of particles in the inner and outer rings is plotted in Fig. |3Jb). The 
average distance between particles in the inner ring is 1.48 ± 0.02 mm, and 1.60 ± 0.03 mm 
for the outer ring. An m = 2 modulation in the particle spacing is visible. 

The particle positions undergo small amplitude thermal fluctuations, mostly due to their 
interaction with the neutral gas component [TS]. For the 1-ring [Fig. [2(a)], velocity distribu- 
tions in the longitudinal (azimuthal) and transverse (radial) directions were computed from 
the time histories of the measured particle positions [23] . The velocity distribution functions 
indicate a minimum measurable velocity of ~ 0.06 mm/s, where a velocity of 1 pixel/frame 
corresponds to v — 0.238 mm/s. Except for a small neighborhood around v — 0, the distri- 
butions in both directions are well fitted by a Maxwellian with a temperature T = 410 K, 
somewhat below 530 K reported for a 2D lattice [23J. This temperature indicates that the 



7 



ring is stable and in equilibrium with the neutral gas. 

The dependence of the power spectral density for the thermal fluctuations vs wavenumber 
k for the 1-ring were computed for the longitudinal (acoustic) and transverse (optical) modes. 
First, the average position of each particle was computed by averaging over all 1025 frames. 
Then, for each frame the longitudinal and transverse displacements were computed and 
Fourier analyzed using a discrete Fourier transform (DFT) to give the amplitude and phase 
of each k vs time, which assumes the equilibrium lattice spacing is independent of position. 
The power spectral density (psd) of the time history of each wavenumber was then computed 
using the maximum entropy method (MEM) [21]. The psd's from four data sets were 
averaged to give the final psd. Since each normal mode acts as an independent damped 
harmonic oscillator, the psd peaks near the mode's natural frequency and so the maximum 
of the psd should follow the dispersion relations [IT] . 

The power spectral density for the acoustic modes is displayed in Fig. ffla). Here the 
horizontal axis runs from mode number to 34, or equivalently ka = to n. The vertical 
axis starts at 0.5 Hz because of low- frequency noise and runs to the Nyquist frequency of 3.75 
Hz. The power spectral density for the transverse mode is shown in Fig. [4](b) . Here we see 
a typical optical mode, with a finite frequency for k — 0, which should be coo [Eq. (|2|)], and 
then a decrease in u with k indicating a backwards wave. The m = 1 mode displays several 
narrow band features due to minute mechanical vibrations coupled to the vacuum chamber 
from our belt-drive mechanical pump. In this case, the electrode and camera oscillate while 
the dust ring remains stationary. 

We fitted the psd for each mode with the equation for a driven damped harmonic oscillator 
for both the longitudinal and transverse directions to determine the natural mode frequency 
bji vs k. Here toi may be slightly above the maximum of the psd because of finite damping, 
where the damping rate 7 ~ 2 s _1 . We then fitted both measured dispersion relations 
simultaneously to theory [Eqs. (|8]) and (J9|], where there are only three free parameters: 
uq, a and R. From this procedure, we found u$ = 16.6 rad/s, a = 1.41 and R = 1.32 ± 
0.05. We believe that this is the first report of very good agreement of experimental data 
simultaneously with both branches of the dispersion relation, confirming that this system 
is well described by a model for an unbounded ID chain of particles interacting through a 
Yukawa potential. For R — 1.32 the critical value of the lattice parameter a = a/r^ for the 
zigzag instability is a c = 1.13. For the ring studied here, a = 1.41 is significantly above 
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Figure 4: (Color online) (a) Power spectral density (psd) for thermally-excited longitudinal modes 
vs mode number and frequency and (b) psd for transverse modes. Purple corresponds to the least, 
and red to the greatest, power density. Open circles are the natural frequency for each mode, and 
the dashed line is the fit to the theoretical dispersion relations. 

a c . Dimensional parameters are tq = 0.77 mm, Xd = 0.82 mm and q = —1.24 x 10 4 e. The 
values of Ad and q are consistent with those found in DONUT for similar plasma conditions 

[3 HMD- 

In summary, we have created ID and quasi-lD dusty plasma rings experimentally, verify- 
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ing previous predictions for charged particles confined in a 2D annular well [131 EH]- Measured 
longitudinal and transverse dispersion relations in the 1-ring exhibit excellent agreement with 
dispersion relations in an unbounded Yukawa lattice, confirming that these systems are well 
described by a model of n identical charged particles confined in an annular potential and 
interacting through a Yukawa potential. These rings provide a new experimental system for 
studying static and dynamic properties of ID and quasi-lD systems. 
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